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ABSTRACT 

Prompted by various analyses of long (Type II) GRB rates and their relationship to 
the cosmic star-formation history, metallicity and luminosity function evolution, we 
systematically analyse these effects with a Monte Carlo code. We test various cosmic 
star-formation history models including analytical and empirical models as well as 
those derived from cosmological simulations. We also explore expressions for metal- 
licity enhancement of the GRB rate with redshift, as presented in the literature, and 
discuss improvements to these analytic expressions from the point of view of galactic 
evolution. These are also compared to cosmological simulations on metal enrichment. 
Additionally we explore possible evolutionary effects of the GRB rate and luminosity 
function with redshift. The simulated results are tested with the observed Swift sample 
including the L, z, and peak flux (log A^ — log P) distributions. The observational data 
imply that an increase in the GRB rate is necessary to account for the observations at 
high redshift, although the form of this enhancement is unclear. A rate increase due 
to lower metallicity at higher redshift may not be the singular cause and is subject 
to a variety of uncertainties. Alternatively, evolution of the GRB luminosity function 
break with redshift shows promise as a possible alternative. 
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Since firmly establishing the cosmologi cal nature of Type II 
(long-soft) gamma-ray b ursts (GRBs) (|Metzger et al. 19971 : 
Ivan Paradiis et al. 19971 ). there have been many predictions 
as to how early in cosmic history GRBs are created. Red- 
shifts for GRBs have been detected more effectively since 
the 2004 launch of the Swift satellite l|Gehrels et al. 20041 ) 
which has the advantage of providing prompt localizations. 
This combined with the dedicated work of ground-based as- 
tronomers has shown progress in pushing toward the theo- 
retical detection limit of about a z ~ 20 (|Abel et al. 200 j : 
iBromm fc Loeb 20021). Record-breakin g bursts, such as 



It is believed that Type ifl GRBs are a product of the 
core-collapse of massive stars, stemming from the evidence of 
an a ssociation of these GRBs with c ore-collapse supernovae 
( Stanek et al. 20031 : llTiorth et al.YoOa ). These observations 
lead naturally to the expectation that the rate of these 
objects follow the cosmic star- fo rma tion his tory (SFH) 
llWiiers et al. 19981: iTotani 19991: iLam b & R eichart 200o( : 



iBlain fc Nataraian 2000l : iPorciani fc Madau 200ll ). Varl 



050904 (z 



Kawai et al. 



20061: 



Cusumano et al. 20061 : iHaislip et al. 20061 : 
Frail et al. 20061 ). 



ous studies have shown that the rate of GRBs does not 
strictly follow the SFH but is actually en h anced by some 
mechanism at high-2 
Guctta & Piran 20o3; 
Kistler et al. 20091: I 



08 0913 {z =6.7 

Greiner et al. 20091 ) and 090423 (z=8.2, ISalvaterra 20091 : 
Tanvir et al. 20091 ) demonstrate just how far these objects 



: (iDaigne et al. 2 006': "Le et al. 200J 
iLi 20081: Kistler ct al. 2008 
Salvaterra 20091: ISalvaterra 2009b 



Salvaterra fc Chincarini 2009 



Qin et al. 201C 



can be detected and warrant a discussion on how bursts that 
occur at such drastically different times in the evolution of 
the universe may or may not differ. 



Wanderman fc Piran 2010l: "Icampisi et al. 2010l) . " be it 
metallicity effects ()Langer fc Norman 20061 : iLi 20081 ). 
selection effects or an increase in luminosity. 

In this analysis we combine and expand various ele- 
ments from these works to further analyse possible GRB rate 
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bee IZhang et al. 20091 for a full discussion on the classification 
of GRBs and a full description of the distinction between Type I 
and Type II bursts. 
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enhancements with redshift and the underlying causes and 
forms of these evolutions using the available observational 
data together with a Monte Carlo code. We look into the 
underlying form of the cosmological SFH, including mod- 
els derived from cos mological smooth-particl e hydrodynamic 
(SPH) simulations (|Choi fc Nagamine 2010I ). metallicity ef- 
fects, rate evolution with redshift and evolution of the break 
luminosity of the GRB luminosity function. 

In Section 2 we present and explain the details of the 
various simulations that were conducted, broken down by 
the form of the SFH or high- 2; enhancement (e.g. metallic- 
ity or evolution effect). Section 3 explains the method of 
testing for consistency. Section 4 details the results for the 
simulations in the same framework as Section 2, dedicating 
a section to each form of enhancement. We conclude with a 
summary and discussion in Section 5. 



2 THEORY AND SIMULATIONS 

One of the major goals of this and previous analyses is to 
constrain the intrinsic distribution of GRBs by utilizing the 
available observed data. We develop a Monte Carlo code that 
randomly creates a set of GRBs, defined by a luminosity and 
redshift pair from assumed luminosity and redshift distribu- 
tions, and cycles them through a series of filters that act as a 
'detection'. The setup is similar to that of lVirgili et al. 20091 
but with various additions and improvements tailored to this 
specific problem. The set of generated bursts is then com- 
pared to the current observations of the luminosity, redshift, 
and peak fiux (i.e. log A'' — log P) distributions. The observed 
GRB rate follows the form 



dN 



Rgrb{z) dV{z) 
dz 



(1) 



dtdzdL 1 + z 

where the {1 + z) factor accounts for the cosmological time 
dilation, Rgrb(z) is the GRB volume event rate (in units 
of Gpc"'' yr~^) as a function of z, $(!/) is the luminosity 
function, and dV{z)/dz the comoving volume element given 
by 



dV{z) 
dz 



Ho {1 + z)^[^Im{1 + z)3 + Qa]^/^' 



(2) 



for a flat A cold dark matter (ACDM) universe. We assume 
Ho — 71 km s~^ Mpc~^, f2m=0.3, and f2A=0.7 throughout. 

Equation [1] has two unknowns to be explored, the 
luminosity function term being the more straightforward. 



texts tSchmit et al. 200ll: Llovd-Ronning et al. 20021: 


Norris 20021: IStern et al. 20021: 


Guetta et al. 20051: 


Llovd-Ronning et al. 2004 


Daigne et al. 20061: 


Coward 20051: ICobb et al. 20061: 


iPian et al. 20061: 


Sodcrbcrg ct al. 20061: IChaoman et al. 20071: 


Liang et al. 2007: iDai 20091: 


fVirgili et al. 20091: 


Qin et al. 2010l: IWanderman & Piran 2010h. We adopt 



the generally accepted broken power-law model: 



(3) 



where a\ and a2 are the power law indices. Lb the break 
luminosity, and $0 a normalization constant. We consider 
solely 'classical' high-luminosity GRBs, ignoring local low- 
luminosity events and their contribution to the luminos- 



ity function jCoward 20051 : iLiang et al. 20071 : iLe et al. 20071 : 
IVirgili et al. 2009 ) in order to have as unbiased a sample as 
possible. 

The GRB rate, Rgrb{z), is the main focus of this analy- 
sis as it is a convolution of the star-formation history, metal- 
licity and evolution effects. Next, we present the specifics of 
each SFH model and enhancement. 

2.1 Cosmic star- formation history 

The cosmic star-formation history is the basis for the rate 
distribution from which we choose our redshift values for 
the simulated GRBs. Many forms are available in the lit- 
erature but it is generally believed that the SFH increases 
rapidly to about z ^ 1 — 2 then slowly falls off toward higher 
redshift, and we use a variety of forms presented in the liter- 
ature. Hopkins and Beacom (2006, or 'HB') have compiled a 
widely acce pted model fit from nume rous multi-band obser- 
vations (see lHopkins fc Beacom 20061 and reference therein). 
Bromm and Loeb (2006 or 'BL') present a comprehensive 
model for the SFH based on a flat ACDM cosmological 
model with the added contribution of Population HI stars 
at high redshift. As a control we also include the SF2 model 
of Porciani and Madau (2001, or 'PM'; See Figure [T| based 
on estimates from UV-continuum and Ho emission. A list 
of SFH models used are summarised in Table 1. 

In addition, we utilize a model derived from cosmo- 
logical SPH simulations of Choi & Nagamine (2010, or 
'CN'). They developed a modified version of GADGET-3 
code (originally described in lSpring cl 2005 ), including radia - 
tive cooling by H, He, and metals I Choi fc Nagamine 20091 ). 
heating by a uniform UV backgroun d of a modified 
Haart and Ma dau (1996) spectrum (iKatz et al. 199(^ : 
iDave et al. 1999"), a sub- resolution model of multiphase ISM 
(Sprin gcl fc Hcr nquist 2003), the "Pr essure" star f ormation 
model (|Schave et al. 2010l : IChoi fc Nagamine 2010l ). and the 
"Multicomponent Variab le Velocity" galactic wind model 
(|Choi fc Nagamine 201lh . They have shown that the metal 
line cooling enhances star f ormation across all redshifts 
by about 10 - 30% (|Choi fc Nagamine 2009), and that 
the Pressure SF model suppresses star formation at high- 
redshift due to a higher t hreshold density for star-formation 
(IChoi fc Nagamine 2010l) with respect to the earlier m odel 
by ISpringel fc Hernguist 200a . IChoi fc Nagamine 201 ll also 
showed that the MVV wind model, which is based on 
both momentum-driven and energy-driven galactic winds, 
makes the faint-end slope of the GSMF slightly shallower 
compared to the constant velocity galactic wind model of 
ISpringel fc Hernguist 20()3l . The adopted cosmological pa- 
rameters are consiste nt with the WMAP best-fit values 
l|Komatsu et al. 201 ll ): Q.m = 0.26, = 0.74, 0.^ = 0.044, 
h = 0.72, Us = 0.96, and erg = 0.80. These values are 
only slightly different from those presented in § 2. The 
results from three simulations with box sizes of comov- 
ing 10,34 and lOO/i"^ Mpc were combined to obtain a 
full SFH, including galax y stellar masses above 10^ Mq 
(|Choi fc Nagamine 20l"Tbl ). 



2.2 Metallicity 

One popular explanation for possible enhancements of 
the GRB rate is the effect due to decreasing metallicity 
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Table 1. Summary of star formation- history models 



SFH Model 


Reference 


PM 


Porciani & Madau (2001) 


HB 


Hopkins & Beacom (2006) 


BL 


Bromm and Loeb (2006) 


CN 


Choi & Nagamine (2010) 



with redshift jFvnbo et al. 20031 : IConselice et al. 20051 : 
Gorosabel et al. 20051 : IChen et al. 2005| : 

Starling et al. 20051: iLanger fc Norman 20061: iLi 200S " 



Qin et al. 2010l : iButler et al. 2010l : ICampisi et al. 201ol ). 

If GRBs occur more frequently in low-metallicity en- 
vironments, then this could be a possible mechanism 
for enhancing the GRB rate at high redshift. Langer & 
Norman (2006; LN) proposed an analytical form for the 
meiss density fraction in galaxies with a mass less than 
M, based on the galaxy stellar mass function (GSMF). 
This mass is then related to the a mount of metals via 
the galaxy mass-me tallicity relation (|Tremonti et al. 20041 : 
ISavaglio et al. 20051 ). givmg: 



^,Z_\^ f[aG + 2,(Z/^o)^10°-^^'^1 
Ze J r[QG + 2] 



(4) 



where f and F are the incomplete and complete gamma 
functions, /3 the power-index of the galaxy mass-metallicity 
relation, and aa the faint-end slope of the GSMF. To be- 
gin, we utilize this relation with constant parameters from 



the literature (ac 



- 1-16. 



-2, e = (Z/Zq) = 0.1; 



iLanger fc Norman 20061 : iLi 20081 ) and then build upon it to 
formulate a weighted version to account for the more realis- 
tic case of variations in metallicity from e = 0.1 — 0.4 (Figure 
[2|. This function also contains many assumptions about the 
underlying mass distribution and mass-metallicity relation 
which we will expand on in §4.4. 

The cosmological simulations derive star formation 
rates for populations of stars from various metallicities with- 
out the need of an external expression. At every time step, 
star particles are created in high-density regions that exceed 
the threshold density according to the star-formation law 
matched to the locally observed Kennicutt (1998) law. Once 
a star particle is created, instantaneous recycling is assumed, 
and the metals are ejected with an yield of Y = 0.02 and 
distributed to the nearby environment by a galactic wind. 
Niino et al. (2011) have used similar simulations to exam- 
ine the metallicity of GRB host galaxies, and found good 
agreement with observations. See Figure [2] for a comparison 
of Equation [4] and the result from cosmological simulations. 

2.3 Rate evolution 

Motivated by the literature (e.g. iKistler et al. 200^ : 
iKistler et al. 20091 : iQin et al. 2010l ) we include a discussion 
on an increase in the GRB rate with redshift as {1 + z)^ . 
This is not to be confused with the evolution of the break 
luminosity of the luminosity function, which has the same 
functional form (See §2.4). The latter has a more physical 
meaning (e.g. GRBs becoming brighter with increased red- 
shift), whereas the former is more of a general statement 
of the GRB rate, increasing in this fashion due to an un- 



specified process. The former may be related to an evolving 
stellar initial mass function with redshift that causes a shift 
to a top-heavy stellar IMF (I Wang et al. 201 ll ). These simu- 
lations assume Rgrb oc SFH x (1-1-2)'' with a non-evolving 
luminosity function. The results of these simulations are pre- 
sented in §4.5. 



2.4 Lt evolution 

The last form of enhancement of the GRB rate is evolution of 
the break of luminosity function. Lb, with redshift. We take 
a similar functional form to the rate evolution, assuming 
that the luminosity function break evolves as (1 -I- z)'^ . This 
increase manifests itself as an increase of bright bursts at 
higher redshifts, which increases their detection rate. Unlike 
the previous section, these simulations assume that Rgrb 
is proportional to the SFH and that the luminosity function 
break, Lt, evolves as (l+z)'' . The results of these simulations 
are presented in §4.6 

2.5 Threshold and other details 

Once a luminosity and redshift pair is chosen according to 
the distributions discussed above, it is necessary to adopt 
a threshold condition that mimics the detector in question. 
We adopt the threshold condition based on the probability 
of triggering Swift derived by Qin et al. (2010): 



5.0P■^■**^ 

0.67(1.0 - 0.40/P)'^■^^ 



P < 0.45 
P > 0.45 



(5) 



where P is the photon flux of the burst in the 15—150 keV 
band. This equation is based on the similarities between the 
BATSE and Swift photon flux samples. By comparing the 
relative number of bursts, both triggered and un-triggered, 
that occur in a particular photon flux bin to the total num- 
ber of bursts it is possible to obtain a probability for trig- 
gering that instrument. Qin et al. (2010) were able to flt the 
distribution to derive Equation [5] When comparing different 
detectors it is necessary to have a normalized value for P, 
and Qin et al. found that the normalization for both detec- 
tors is similar, so one expression is a good approximation for 
both detectors. For comparisons with the observed redshift 
and luminosity samples, we include an additional probabil- 
ity for the detection of a redshift since not all bursts have 
redshifts. Similarly, a probability of assigning a redshift is 
found by looking at the distribution of bursts with redshift 
versus the total number of bursts per redshift bin. Qin et 
al. (2010) did not find a large distinction between these two 
samples but parametrize the probability as 

1.61 logP 



ri^ = 0.26 + 0.032e" 



(6) 



Both expressions and the log A'^ — log P analysis depend on 
the calculation of the photon flux. The energy flux is cal- 
culated directly from Fpk — L/4,-KD\{z)k where Dl{z) is 
the luminosity distance at a given redshift, and k is the k- 
correction 



k = 



/(1+^) 



EN{E)dE 



m EN{E)dE 



(7) 



which corrects the flux from the bolometric 1 — 10* keV 
band into the observed band (ei,e2). In Equation [71 N(E) 
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is the photon spectr um of the GRB, which we assume to 
be a Band function l|Band et al. 199a ). The Band function 
is a smoothly joined power law function that has pre- and 
post-break slopes a and /? around a break energy Eo- The 
peak of the fF,^ spectrum, Ep^ak, is related to this energy 
by Epeak = {2 + a)Eo. Since the peak energy of bursts 
change s with the energy of the burst (i.e. the Amati re- 
lation; lAmati et al. 20021 : iLiang fc Dai 20041 ') we utilize the 
relation derived by Liang et al. (2004) to assign values for 
the simulated Ep^ak- 

Speafc/200keV = C-'(L/10'"erg s.-^f'^, (8) 

where C is a random uniform deviate between [0.1,1] and L 
the luminosity. We also randomly sample a and /3 uniformly 
between -0.83 < a < -1.2 and -2.1 < P < -2.5, which 
roughly correspond to the observed limits of these values. 
From these spectral parameters and peak energy flux, F, we 
then calculate the peak photon flux in the detector energy 
band (ei, 62) via 

F N(E)dE 

p EN{E)dE ■ ^ ' 



3 TESTING FOR CONSISTENCY 

Once we have a set of simulated bursts that are 
'detected' and follow the GRB rate and luminos- 
ity function described above, it is necessary to test 
the consistency with the observed data. Our sam- 
ple consists of 166 Swift and HETE era GRBs with 
known redshift thro ugh September 20 9. We remove the 
Type I GRBs ( see 'Zhan g et al. 20091 : IVirgih et al. 201ll : 
iKann et al. 201 ll ). outlying low-luminosity bursts (GRB 
980425, 060218), and any bursts with disputed or non-secure 
redshifts. This assumes that the observed redshift sample 
is the true intrinsic sample, although the detection of red- 
shifts depends on a variety o f observation al factors and po- 
tential biases (|Fiore et al. 2007 : Jakobsson et al. 201 ll ). In 
order to calculate the bolometric peak luminosity it is nec- 
essary to have the energy or photon flux as well as the 
spectral information for the fc-correction. Most Swift bursts 
are flt by a simple power-la w spectrum due to sm all band 
pass of the Swift detector (|Sakamoto et al. 2007^ . We in- 
stead assume all bursts to have a Band function spectrum 
(|Band et al. 19931 ) utilizing the observed pre-break slope 
and assuming the typical value of 2.5 for the post-break 
slope, as assuming a simple power-law extending to high en- 
ergies will inevitably overestimate the high-energy contribu- 
tion to the flux. This assumption is validated by observations 
of GRBs observed jointly with Fermi LAT and GBM, where 
a Band function spectrum is seen over m any orders of mag- 
nitude and exten ding to high energies (|A.bdo et al. 20091 : 
IZhang et al. 20l"ll ). Epeak is used from the literature unless 
absent, in which case the value from catalogue of Butler et 
al. (2010) is used. 

Next, the simulated set of 175 simulated bursts is 
compared to the observational sample with the k-sample 
Anderson-Darl ing (AD) test for consis tency between two 
distributions (|Scholz &: Stephens 19861 ') ■ giving the first 
three constraints (L alone, 2 alone, L ~ z together) to our 



models. Each criterion gives a contour showing the consis- 
tency with the observed sample in the (ofi, Lt) plane (Fig- 
ure 3^) with 0:2 a constant of 2.2 or 2.5. We test both 
values of 0:2 and results indicate which slope was used. 
Results are generally insensitive to the choice of a 2 and 
we take values quoted in the lite rature (|Liang et al. 20071 : 
IVirgih et al. 20091 : lOin et al. 2010l ). 

The two logA'^ — logP tests are also conducted 
with the AD test and compares the simulations to the 
CG-RO/BATSE and Swtft photon flux samples. In order to 
have the most unbiased and complete sample from each, 
they are truncated at 0.4 [50-30 ke'V] and 1 ph cm~'^ s~^ 
[15—150 ke V] , respectively (see iLaredo fc Was erman 19981 : 
iBand 20061 ') and compared to the 1143 BATSE and 380 trig- 
gered Swift bursts. A summary of various models and sta- 
tistical results are presented in Tables 2-4. 

4 RESULTS 

4.1 GRB rate oc SFH 

Clearly the simplest scenario possible for the GRB rate, this 
set of simulations showed little consistency with the current 
observations. Out of the four possible SFH models (Porciani 
& Madau (2001), Hopkins and Beacom (2006), Bromm and 
Loeb (2006), Choi and Nagamine (2010)), only the Bromm 
and Loeb model showed consistency with the observations 
(Table 2). The luminosity function is generally constrained 
to be shallow, with pre-break slopes generally more shallow 
than ~ —0.2, and shows that there is the need for some 
form of increase of the rate compared to other SFH mod- 
els at higher redshifts. Determining the form and possible 
cause(s) of this increase in GRB rate is a major goal of this 
analysis. This model, however, is based on the theory that 
the rate enhancements produced at high z are attributed 
to the contribution of Population HI stars that were devel- 
oping in the early universe z ~ 15. To date, observations 
have not shown that GRBs a rise from Population HI stars 
(ITanvir et al. 20091 : ISalvaterra 2009. ') and we caution draw- 
ing an association based solely on the form of this SFH. 



4.2 GRB rate oc SFH+Metallicity cutoff 

The next step, as detailed above, is to consider the addition 
of a term that in some manner accounts for an increase in 
the GRB rate in lower metallicity environments. A key com- 
ponent believed to aid in the creation of Type II GRB is a 
fast-rotating core of the progenitor star. Low metallicities 
may help in reducing the mass-loss rate and retain sufflcient 
angular momentum to keep the star rotating quickly and 
assist in the formation of the GRB jet. First we con sider the 
formalism of L anger and Norman (2006; see also iLi 20081 : 
iQin et al. 20101 ) as detailed in Eqn |4j for all models with 
the exception of CN, as they form stars self-consistently ac- 
cording to metal line cooling rates in the simulation. The 
derivation is straightforward and clear in Langer and Nor- 
man (2006), however, there are various assumptions in this 
model that need to be addressed. The basis for this rela- 
tion is the galaxy stellar mass function (GSMF), which is 
assumed to be a Schechter function (i.e. a power law with 
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Figure 1. Panel (a): Star-formation history models utilized in 
the analysis. Panel (b): Simulated relative number of GRBs per 
unit co-moving volumex (l-|-z)^^ for different models. This panel 
shows the output from the code without a threshold, so as to 
check the underlying distribution and see the relative affects of 
the metallicity relations on the base SFH. 

exponential cutoff; Eqn 1 in lLanger fc Norman 20061 and ref- 
erences within). The amount of galaxy stellar mass within 
a mass M is then related to the amount of metals by the 
mass-metallicity relation, of the form MjM* = K{Z/Zq)'^ , 
where K and P are constants that are constrained by obser- 
vation l|Tremonti et al. 20041 : ISavaglio et al. 20051 ) and M* 
is the characteristic mass of the GSMF. Previous stud- 
ies do not address the scatter and/or evolution with red- 
shift of the GSMF faint-end slope ac, and assume that 
the average cosmic metallicity of the Universe decreases 
as d[Z]/dz = —0.15 dex. The normalization of the mass- 
metallicity relation changes with redshift and for a particu- 
lar metal cut e = (Z/Zq), the mass fraction of metals also 
changes. This effect is reflected in the IQ''-^^''^ term of Equa- 
tion 3] Changes to this term are not examined explicitly. 

Using this relation with the parameters assumed in 
Langer and Norman (2006) (ac ~ —1.16, P — 2) and a cut- 
off metallicity for production of GRBs of e = Z/Zq = 0.1, 



Figure 2. Fractional stellar mass density contained in galaxies 
with metallicities below Z/Zq (Eq.|4]l from Langer and Norman 
(2006, LN), including different metal cuts {Z/Zq = 0.1,0.4) and 
modifications from weighting. The expression derived from the 
star-formation history of CN with a metal cut of 0.1 is included 
for comparison. 



we find that no models agree with the L and z constraints to 
the 2a level. The cosmological simulation results are similar 
and show that a strict metal cut at Z/Zq = 0.1 is insufficient 
to explain the observations. 



4.3 GRB rate oc SFH+Weighted metallicity 

Building upon the previous section, we introduce a weight- 
ing to the metallicity cut in order to broaden the scope of 
Eqn 131 It is more realistic to consider a range of metallicities 
in which GRBs can occur, especially since GRBs have been 
observed in environments with metallicity greater than 0.1 
(|Holland et al. 2010l : iLevesaue et al. 20101 ). Instead of tak- 
ing the value of the relation from Langer and Norman (2006) 
for a particular value of e, we instead weight the effect of dif- 
ferent metallicities, ranging from Z/Zq = 0.1 — 0.4, for a par- 
ticular redshift. We have utilized a Gaussian (with cr = 0.1 
or 0.2) to weight the contributions of metals so that there 
is an exponential (rather than sharp) cut above the critical 
metallicity. The contributions from various metallicities are 
then added with proper weighting to produce an 'effective' 
^' (Equation |3| . This approach yields an intermediate so- 
lution between strict metal cuts of 0.1 and 0.4 (Figure [2]). 
A similar approach is taken for the CN model utilizing the 
SFR for different metal cuts provided from the simulation 
instead of applying Eqn|4l 

Using this formulation we re-run the previous SFH mod- 
els and find that the HB model is the only model that can 
pass all of the observational tests, including the log A'^ — 
log P, giving luminosity function parameters in the range of 
(ai, 02 = 2.2, Lb) = (0.11 - 0.19, 2.2, 6 - 10 x 10^^ erg s~^) 
(Figure 21 Table 2). The BL model, with its intrinsically 
large rate at high-z, overproduces bursts at high-z when 
metallicity is added. 
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4.4 More on the metallicity approximation 



Up to this point, the analysis does not directly compare 
how the assumed metallicity relation (Eqn |3| affects the 
GRB rate compared to the cosmological simulations. This 
is an important and related topic, since the models that 
use the LN expression show consistency with the obser- 
vations of Type II GRBs, while the more rigorous and 
complete method of CN to calculate the metallicity shows 
no consistency. The differences, we come to find, are non- 
negligible. Why would the HB and CN models, whose total 
star-formation rates are quite similar, differ so largely when 
their respective metal cuts are applied (Figure [T]d)? The re- 
lation from Langer and Norman (2006) is an approximation 
to a very complex problem in galaxy evolution. The cosmo- 
logical simulations by Choi and Nagamine (2010) address a 
variety of effects that contribute to the metal distribution 
(such as mixing due to galactic outflow and tidal disrup- 
tion), and calculate the star-formation rate selfconsistently 
according to varying metal line cooling rate. From those val- 
ues a realistic view of how the total rate is affected by the 
reduction in metallicity can be calculated, which is just what 
Eqn [4] shows: the net effect to the total star-formation rate 
by a metallicity cut at Z/Zq. The curves for various values 
of Z/Zq are shown together with the equivalent expression 
from the CN (Figure [2] and [S]). 

Since these expressions are so different, we attempt to 
look at the structure of the Langer and Norman (2006) ex- 
pression and see if any part(s) can be improved to create 
a more realistic view of how metallicity affects the rate of 
GRBs. The first major assumption in Eq. [4] is the constant 
value of the G SMF slope, ag, which is observed to be steep- 
ening with z l|Bouwens et al. 2010l and references therein), 
suggesting a larger number of lower-luminosity galaxies at 
higher redshifts. 

Bouwens et al. (2010) detailed several observations of 
galaxies at 2 ~ 7 — 8 and summarised the evolution of the 
luminosity function of galaxies. From their Figure 15 we are 
able to extract the slope of the GSMF as a function of red- 
shift to incorporate into our code. Using a spline fit and 
cubic interpolation we are able to approximate the behav- 
ior of qg both at the data point as well as maximum and 
minimum values from the error bars provided, which range 
from about — 1 > ag > —2 in the range 2 ~ — 10 (Figure 
[5}. We consider only values above ag = — 2 as the metal- 
licity relation is undefined at the value ac 4-2 = 0, which 
affects the minimum error bar approximation. For that case 
we assume ac is constant with the value of the lowest data 
point (-1.99) out at higher redshifts. Above 2 = 8 we again 
assume all values of ac are constant. As shown in Figure 
[31 the evolution of ac implies a faster cosmic metallicity 
enrichment than just applying the unaltered expression and 
pushes the curve toward lower redshift and closer to the re- 
sults from cosmological simulations. Using the values of the 
upper error bars of Figure[3^ gives solutions that are similar 
to the weighted expression of LN (Table 3). 

This more realistic approach show some consistency 
with observation, but only in a few cases. The HB model uti- 
lizing the upper limits of the evolution of ac and a metal cut 
of Z/Zq =0.1 shows some consistency with the L and z con- 
straints and consistency with all logA'^ — logP constraints, 
while models with metallicity weighting show results only 
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Figure 3. Panel (a): Redshift evolution of the GSMF faint-end 
slope, Qfj, including error bars (Bouwens et al. 2010). Panel (b:) 
Effect of the evolution of ag on the expression from LN. The 
Z/Zq = 0.1 cut expression from LN and CN are included for 
comparison. 



in the 3cr contour. All other HB models show little consis- 
tency in all tests while the BL model shows some areas in 
the L and z constraints but shows a large deviation in the 
log iV - log P resuhs for both BATSE and BAT. The lack of 
consistency that is evident in most models is generally at- 
tributed to an overproduction of bursts at z ~ 1 — 2. Models 
utilizing the original framework of LN may show more con- 
sistency with the observations, but these set of updates are 
a promising and necessary direction for study that, with fur- 
ther enhancements, might be able to fully explain the rate 
increase. 



4.5 GRB rate cx SFH X (1-1-2)* 

As detailed in the literature (iKistler et al. 20081 : 
iKistler et al. 20091 : lOin et al. 20101 ) we also consider 
an increase in GRB rate as (1 -I- 2)'', where 3 — 0.2, 0.5, 0.8. 
We consider all of the SFH models with and without 
metallicity enhancements (no GSMF evolution) and find 
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that a few of these models are able to pass the L and z con- 
straints but fail to pass the BATSE and Swift log A'^ — log P 
constraints (Table 4). 

4.6 Luminosity function break evolution oc (1 + z)'' 

Lastly we consider the evolution of the luminosity func- 
tion break as detailed above. We see some consistency with 
the CN model and evolution with 7 ~ 0.5 — 1.5. The 
3a regions for the CN models show areas of consistency, 
with a few showing 2cr significance (i.e. 7 = 1.0, 1.3) (Fig- 
ure [5|. The general trend is again for shallow luminosity 
function slopes, the best models occurring in the area of 
(01,02,1/6,7) = (0.5,2.2,3 X 10^^ erg s~\ 1.0). The HB 
models show some consistency to So" in the same regions, 
although not as broadly as the CN model (Table 4). 



5 SUMMARY AND DISCUSSION 

Our work supports the idea that the GRB r ate is enhanced 
at higher redshift (iDaigne et al. 2009: iLe et al. 2007 



Li 20081: iKistler etld 



[Salvate 
iQinc 



2008 



ISalva terra 2009b 
ct al. 201c 



Guetta fc Piran 20071 : 

Kistler et al. 20091 : [ Salvaterra 20091 ; 

Salvatcrra & Chi ncarini 20091 : 

Wanderman fc Piran 2010l '). The form of this increase, 
however, is still unclear. We have tested various SFH 
models and enhancements to the GRB rate, reflecting 
possible effects from changing cosmic metallicity and other 
evolutionary effects, with a Monte Carlo code. The resulting 
output was then tested for consistency with a variety of 
available Swift and BATSE data, including the L, z, and 
peak photon flux distributions. Even when considering 
a numerical simulation model that takes into account a 
variety of realistic galactic evolution effects, both with 
and without me tal cuts, and a metallici ty relation based 
on the GSMF (jLanger fc Norman 20061 ) our models do 
not show strong consistency with the observed sample, 
although we believe this is the right direction for this type 
of study. This may indicate that metallicity is not solely 
responsible for the increased rate and that perhaps some 
other type of enhancement is needed. To this end, we test 
both GRB rate evolution and luminosity function (break 
luminosity) evolution with redshift, finding that the latter 
is allowed within the constraints of the BATSE and Swift 
data with moderate (oc Lf, x (1 -f z)^'^'^~^''^) evolution. This 
statement has, of course, a few caveats. Embedded in the 
metallicity relation are a variety of assumptions about the 
GSMF and the observed mass-metallicity relation. Laskar 
et al. (2011) show, using HST observations of GRB host 
galaxies, that the metallicity relationship likely evolves 
between redshifts of 3—5, which would further affect the 
results. It is possible that other combinations of parameters 
or assumptions might yield a more realistic relation, and 
we suggest further work on how the GSMF and stellar IMF 
work in tandem to affect the problem at hand. In addition, 
recent works have studied the M — Z relation of Type II 
GRBs and found that the hosts lie b elow the SPSS M - Z 
relation (iKocevski fc West 2oTll : iMannucci et al. 201ll : 
ICampisi et al. 201 il l. This adds further evidence to the fact 
that the assumption of this relation for these types of bursts 
is likely not valid, and perhaps a consequence of the active 



star-formation environment instead of a strict metallicity 
cut. We have explored some basic changes, such as the 
evolution of the GSMF faint-end slope, but a comprehensive 
study of this relation or a realistic alternative, are needed. 

We have detailed a numerical and statistical ap- 
proach aimed at understanding the properties of the GRB 
rate in the context of the cosmic star-formation history, 
including the constraints from newly discovered high-2: 
bursts and the possibly effects of metallicity and vari- 
ous types of evoluti on. Recent work s have addressed this 
prob lem in similar (iQin et al. 2010l ) and fully analytical 
( Wa nderman fc Piran 20101 ) ways, and share some common 
points, although both call on fairly strong evolution of the 
GRB rate ((l + .z)* ~ 0.6 — 2) which we do not find. Our work 
also benefits from the i nclusion of a fully nume rical star- 
formation history model (|Choi fc Nagamine 2010l l as well as 
a probing of the metallicity relation and cosmological con- 
siderations that may affect the GRB rate which are not in- 
cluded in contemporary works on the subject. Butler et al. 
(2010) do not find evidence of strong luminosity function or 
GRB rate evolution and find that a smoothed metallicity cut 
of Z/Zq — 0.2 — 0.5, following the metallicity considerations 
of LN, can account for the observations of the current Swift 
sample, although they acknowledge that there are large er- 
rors bars. They also do not include evolution of the GSMF, 
which may account for the differences with this work. In 
addition, we analyse most components separately, and it is 
possible that the observed distribution is a superposition of 
a variety of effects. With enough computational time the 
various combinations of effects can and should be tested. 

By fitting the redshift distribution and log A''— log P dis- 
tribution of BATSE and Swift bursts, Campisi et al. (2010) 
have reached the similar conclusion that Type II GRBs are 
unbiased tracers of the star-formation history. Their anal- 
ysis supports two possible scenarios: (i) a model with no 
metal cuts and a strongly evolving luminosity function or 
(ii) a non-evolving luminosity function with a metal cut of 
Z/Zq < 0.3. Both scenarios assume and fit a Schechter lumi- 
nosity function. This results are similar to the results pre- 
sented here, although the luminosity evolution is stronger 
for the non-metal cut case and the authors claim such large 
changes in GRB properties with redshift as unrealistic, fa- 
voring a model with a metal cut and no luminosity function 
evolution. 



ACKNOWLED GMENTS 

This work is supported by NSF through grant AST- 
0908362, and by NASA through grants NNX10AD48G 
and NNX10AP53G. KN is supported in part by the 
NSF grant AST-0807491, NASA grant HST-AR-12143-01- 
A, National Aeronautics and Space Administration under 
Grant/Cooperative Agreement No. NNX08AE57A issued by 
the Nevada NASA EPSCoR program, and the President's 
Infrastructure Award from UNLV. This research is also sup- 
ported by the NSF through the TeraGrid resources provided 
by the Texas Advanced Computing Center. 



8 F. J. Virgin, B. Zhang, K. Nagamine and J.-H. Choi 



REFERENCES 

Abdo A. A., et al. 2009, Science, 323, 1688. 
Abel T., Bryan G. L., Norman M. L., 2002, Science, 295, 
92 

Amati L., et al., 2002, A&A, 390, 81 
Band D. L., 2006, ApJ, 644, 378 
Band D., et al. 1993, ApJ, 413, 281 
Blain A. W., Natarajan P. 2000, MNRAS, 312, L35 
Bouwens R. J., et al. 2010, arXiv: 100674360 (update?) 
Bromm V., Loeb A., 2002, ApJ, 575, 111 
Bromm, V., Loeb, A. 2006, ApJ, 642, 382 
Butler N., et al., 2010, ApJ, 711, 495 
Campisi M. A., Li L.-X., Jakobsson R 2010, MNRAS, 407, 
1972 

Campisi M. A., Tapparello C, Salvaterra R., Mannucci F., 
Colpi M., MNRAS, submitted fpreprint: larXiv:1105.l"378ll 
Chapman R., et al., 2007, MNRAS, 382, 21 
Chen, H.-W., Prochaska, J. X., Bloom, J. S., Thompson, 
L B. 2005, ApJ, 634, L25 
Choi J.-H., Nagamine, K. 2009, MNRAS, 393, 1595 
Choi J.-H., Nagamine K. 2010, MNRAS, 407, 1464 
Choi J.-H., Nagamine K. 2011, MNRAS, 410, 2579 
Choi J.-H., Nagamine K. 2011b, submitted 

(|arXiv:1101.5656|l 
Cobb B. E. et al., 2006, ApJ, 645, L113 
Conselice, C. J., et al. 2005, ApJ, 633, 29 
Coward D. M., 2005, MNRAS, 360, 77 
Cusumano G., et al., 2006, Nature, 440, 164 
Dai X. Y., 2009, ApJ, 697, 68 

Daigne F., Rossi E., Mochkovitch R. 2006, MNRAS, 372, 
1034 

Dave, R., Hernquist, L., Katz, N., Weinberg, D. H. 1999, 
ApJ, 511, 521 

Fiore F., Guetta D., Piranomonte S., D'Elia V., Antonelli 

L. A.. 2007, A&A, 470, 515 
Frail D. A., et al. 2006, ApJ, 646, L99 
Fynbo J. P. U., et al. 2006, A&A, 451, 3 
Fynbo, J. P. U., et al. 2003, A&A, 406, L63 
Gehrels N. et al. 2004, ApJ, 611, 1005 
Gorosabel, J., et al. 2005, A&A, 444, 711 
Greiner J. et al. 2009, ApJ, 693, 1610 
Guetta D., Piran T., 2007, JCAP, 7, 3 
Guetta D., Piran T., Waxman E., 2005, ApJ, 619, 412 
Haardt, F., Madau, P. 1996, ApJ, 461, 20 
Haislip J. B., et al. 2006, Nature 440, 181 
Holland S. T., et al. 2010, ApJ, 717, 223 
Hjorth J., et al. 2003, Nature, 423, 847 
Hopkins, A. M., Beacom, J. F. 2006, ApJ, 651, 142 
Jakobsson P., Malesani D., Hjorth J., Fynbo J. P. U., 

Milvang- Jensen B. 2011, AN, 332, 276 
Katz, N., Weinberg, D. H., Hernquist, L. 1996, ApJS, 105, 

19 

Kann A., et al. ApJ in press. arXiv:080 4.1959] 

Kawai N. et al. 2006, Nature 440, 183 

Kennicutt, R. C. 1998, ARA&A, 36, 189 

Kistler M. D., Yiiksel H., Beacom J. F., Hopkins, A. M., 

Wyithe, J., Stuart B. 2009, ApJ, 705, 104 
Kistler M. D., Yiiksel H., Beacom J. F., Stanek K. Z., 2008, 

ApJ, 673, 119 
Kocevski D., West A. A. 2011, ApJ, 735, L8 
Komatsu, E., et al. 2011, ApJS, 192, 18 



Lamb D. Q., Reichart D. E. 2000, ApJ, 536, 1 
Langer N., Norman C. A., 2006, ApJ, 638, L66 
Laredo T. J., Wasserman I. M., 1998, ApJ, 502, 75 
Laskar T., Berger E., Chary R.-R., 2011, preprint: 

arXiv:1102:1019 
Le T., Dermer C. D., 2007, ApJ, 661, 394L 
Levesque E., Kewley L. J., Graham J. F., Fruchter A. S. 

2010, ApJ, 712, 26 
Li L. X. 2008, MNRAS, 388, 1487 

Liang E. W., Zhang B., Virgih F., Dai Z. G., 2007, ApJ, 
662, 1111 

Liang E.-W., Dai Z. G., 2004, ApJ, 606, L29 
Lloyd-Ronning N. M., Dai X., Zhang B., 2004, ApJ, 601, 
371 

Lloyd-Ronning N. M., Fryer C, Ramirez-Ruiz E., 2002, 
ApJ, 574, 554 

Mannucci F., Salvaterra R., Campisi M. A. 2011, MNRAS, 
414, 1263 

Metzger M. R., Djorgovski S. G., Kulkarni S. R., Steidel 
C. C, Adelberger K. L., Frail D. A., Costa E., Frontera 
F., 1997, Nat, 387, 878 

Norris J. P., 2002, ApJ, 579, 386 

Niino Y., Choi J.-H., Masakazu A. R., Nagamine K., Totani 

T., Zhang B. 2011, ApJ, 726, 88 
Pian E. et al., 2006, Nat, 442, 1011 
Porciani C, Madau P. 2001, ApJ, 548, 522 
Qin S. F., Liang E. W., Lu R. J., Wei J. Y., Zhang S. N. 

2010, MNRAS, 406, 558 
Sakamoto T. et al, 2007, ApJ, 669, 1115 
Salvaterra R. et al. 2009, Nat, 461, 1258 
Salvaterra R., et al. 2009, MNRAS, 396, 299 
Salvaterra R., Chincarini G. 2007, ApJ, 656, L49 
Savaglio, S. 2006, NJPh, 8, 195 
Savaglio S., et al. 2005, ApJ, 635, 260 
Soderberg A. M. et al., 2006, Nat, 442, 1014 
Schaye, J., et al. 2010, MNRAS, 402, 1536 
Schmidt M., 2001, ApJ, 552, 36 

Scholz, F. W., Stephens, M. A. 1987, J. Amer. Statistical 

Assoc., Vol. 82, No. 399, 918 
Springel, V. 2005, MNRAS, 364, 1105 
Springel, V., Hernquist, L. 2003, MNRAS, 339, 289 
Stanek K. Z., et al. 2003, ApJ, 591, L17 
Starling, R. L. C, et al. 2005, A&A, 442, L21 
Stern B. E., Atteia J.-L., Hurley K., 2002, ApJ, 578, 304 
Tanvir N. et al. 2009, Nat, 461, 1254 
Totani T. 1999, ApJ, 511, 41 
Tremonti C. A., et al. 2004, ApJ, 613, 898 
van Paradijs J. et al., 1997, Nat, 386, 686 
Virgih F. J., Zhang B., O'Brien P. T., Troja E. 2011, ApJ, 

727, 109 

Virgih F. J., Liang E.-W., Zhang B., 2009, MNRAS, 392, 
91 

Wanderman D., Piran T. 2010, MNRAS, 406, 1944 
Wang F. Y., Dai Z. G., 2011, ApJ, 727, 34 
Wijers R. A. M. et al. 1998, MNRAS, 294, L13 
Zhang, B., et al. 2009, ApJ, 703, 1696 
Zhang B.-B., et al. 2011, ApJ, 730, 141 



Gamma-ray burst rate 9 





Figure 4. HB SFH model with tlie weighted LN expression. Panel (a): Contours for consistency in BOTH L and z. Dark grey = 2a 
contour and light grey 3a- contour. Panel (b): Sample 2D distribution from area of consistency in 2a region, (oi, 02, Lf,) = (0.15, 2.2, 8 X 
10^^ erg s~^). Panels (c) and (d): BATSE and Swift \ogN — logP for same parameters as panel (b). 
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Figure 5. CN SFH model with luminosity break evolution ex. {1 + z)^'^. Panel (a): Contours for consistency in BOTH L and z. Dark 
grey = 2a contour and light grey Scr contour. Panel (b): Sample 2D distribution from area of consistency in 2cr region, (01,021^6,7) = 
(0.05, 2.2, 3 X 10^^ erg s~^, 1.0). Panel (c) and (d): BATSE and Swift log — logP for same parameters as panel (b). 
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Table 2. SFH models and test statistics for a variety of simulations, including a metallicity cut of Z/Zq = 0.1. If consistency is found with 
the L and z samples, indicated by a ' Y' in the second column, then the LF parameters are listed with the outcomes of the log N — log P 
analysis. Later models include the addition of metallicity in the form of the expression from LN, but with a Z/Zq = 0.1 cut as well as 
the weighted expression. An 'N' in the significance column indicates that test fails beyond a 3a level. A similar analysis, with similar 
results, was conducted for models with the luminosity function post-break slope 02 = 2.5. Columns 8 and 10 indicate the significance 
level for the BATSE ad Swift log N — log P tests, respectively. 



Model 


L-z? 


LF parameters 


z 


L 


Significance 


BATSE LNLP 


Sig 


Swift LNLP 


Sig 




Y/N 


(ai, if,, 02) 


Stat, P-value 


Stat, P-value 


z/L 


T stat, P-value 




T Stat, P-value 




GRB rate oc SFH 




















02 =2.2 




















HB 


N 


















BL 


Y 


(0.01,500,2.2) 


-0.63727 , 0.51956 


1.49613, 0.07869 


lCT/2a- 


0.92194, 0.13939 


2ct 


0.18829 0.28273 


2a 




Y 


(0.11,600, 2.2) 


1.16638, 0.10932 


-0.34006, 0.42921 


2a-/lCT 


0.42661, 0.22637 


2a 


0.97869, 0.13178 


2a 




Y 


(0.2,900,2.2) 


1.4281, 0.0842 


-0.33604, 0.428 


2ct/1ct 


3.72049, 0.01036 


3a 


0.95844 0.13444 


2a 


CN 


N 


















PM 


N 


















GRB rate oc SFH + Metallicity 
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N 


















BL+Li 


N 


















CN 0.1 cut 


N 


















HB+Li weighted 


Y 


(0.11,600,2.2) 


0.88134, 0.14509 


-0.08952, 0.35632 


2ct/1ct 


1.973, 0.04921 


3ct 


0.17974, 0.28485 


2a 




Y 


(0.19,1000,2.2) 


0.24896, 0.26793 


-0.80114, 0.5692 


2<t/1o- 


0.3929, 0.23235 


2a 


0.48579, 0.21293 


2a 




Y 


(0.15,800,2.2) 


0.29743, 0.25647 


-0.4227, 0.46411 


2ct/1ct 


1.07391, 0.11987 


2a 


-0.21968, 0.39359 


la 


BL+Li weiglited 


N 



















CN weighted N 



Table 3. SFH models and test statistics for models with evolving GSMF faint-end slope, oq, in the metallicity equation of LN. If 
consistency is found with the L and z samples, indicated by a 'Y' in the second column, then the LF parameters are listed with the 
outcomes of the logN — logP analysis, 'sml scatter' and 'Irg scatter' indicate the amount of scatter in the weighting of the metallicity 
relation, 0.1 and 0.2 respectively, 'central values' and 'upper limits' indicate what set of ac values were used in the analysis, those 
corresponding to the data point value or the upper limits in Figure 3a. An 'N' in the significance column indicates that test fails beyond 
a 3(7 level. All models assume the post-break slope of the luminosity function 02 = 2.2. Columns 8 and 10 indicate the significance level 
for the BATSE ad Swift log AT — log P tests, respectively. 



Model 


L-z? 
Y/N 


LF parameters 
(ai,i(,,a2) 


z 

Stat, P-value 


L 

Stat, P-value 


Significance 
z/L 


BATSE LNLP 
T stat, P-value 


Sig 


Swift LNLP 
T stat, P-value 


Sig 


Models including (xq evolution 




















(data point values) 




















HB+Li+ac ovol 
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(0.1,800,2.2) 


2.85189, 0.02164 


-0.61141, 0.51165 
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0.21195, 0.2769 
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0.84384, 0.15056 


2a 
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(sml scatter) 
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HB+Li+cKQ c\'ol+wcighting 




















(Irg scatter) 
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BL+Li+ac evol 


Y 


(0.41,900,2.2) 


0.19157, 0.28191 


-0.37904, 0.44092 


2ct/1ct 


18.10581, 
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1.49981, 0.0784 


2a 


BL+Li+ac evol+weighting 




















(sml scatter) 
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(0.39,800,2.2) 
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2.73904, 0.02396 
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(Irg scatter) 


Y 
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0.74228, 0.16633 


2a 


Models including ctQ evolution 
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HB+Li+ttG evol 
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-0.17686, 0.38088 
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4.42284, 0.00586 
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0.04375, 0.31989 


2a 
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(sml scatter) N 
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(Irg scatter) N 
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(sml scatter) N 
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Table 4. SFH models and test statistics for models with rate evolution proportional to + and luminosity function break luminosity 
(Lj,) evolution as (1 + z)'^ . For the former we show the results for 5 = 0.2 as an example. The results for 5 = 0.5 and 0.8 are similar, 
showing inconsistency with the BATSE and many Swift logN — logP constraints. An 'N' in the significance column indicates that test 
fails to beyond a 3a level. Columns 8 and 10 indicate the significance level for the BATSE and Swift log N — log P tests, respectively. 



Model 
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LF parameters 
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Significance 
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_ 


_ 


_ 
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BL+LN weighted 


Y 


(0.54,900,2.2) 










N 




N 




Y 


(0.46,700,2.2) 
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- 
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N 
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N 




Y 


(0.4, 500, 2. 2) 


- 


- 


- 


- 


N 


- 


N 




Y 


(0.24,400,2.2) 


- 


- 


- 


- 


N 


- 


N 


HB+LN weighted 


Y 


(0.54,900,2.2) 
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N 
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3ct 
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^0.04,900, J. J; 










N 




3ct 


CN weighted 
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LK break evolution 




















i(, CX i(, * (1 + 2)1" 




















7 =1.0 




















HB 


Y 


(0.15,500) 


2.77368 , 0.02321 


-0.37816, 0.44065 


3ct/ Ict 


8.98113, 0.00015 


N 


2.04678, 0.04582 


3ct 


CN 


Y 


(0.05,300) 


0.9696 0.13296 


-0.66268 0.5273 


2ct/ lo" 


0.13296 0.29664 


2ct 


-0.00699 0.33353 


Ict 




Y 


(0.09,400) 


1.33265, 0.0926 


-0.49987, 0.47758 


2o-/ lo- 


2.86931, 0.0213 


3ct 


0.64614, 0.18262 


2ct 




Y 


(0.23,500) 


2.99502, 0.01906 


0.24022, 0.27003 


3o-/ iCT 


-0.59966, 0.50807 


Ict 


0.50419, 0.20924 


2ct 




Y 


(0.23,800) 


1.94944, 0.05035 


1.81983, 0.05714 


2o-/ 2ct 


9.08464, 1.40E-04 


N 


1.49638, 0.07867 


2ct 




Y 


(0.16,600) 


0.97282, 0.13253 


-0.03763, 0.34191 


2ct/ Ict 


6.35583, 0.00124 


N 


2.12331, 0.04256 


3ct 


7 =1.1 




















CN 


Y 


(0.13,300) 


1.1944, 0.1063 


-0.60168, 0.50868 


2ct/ Ict 


-0.10903, 0.36181 


Ict 


0.31864, 0.25158 


2ct 


7 =1.2 




















HB 


Y 


(0.12,300) 


3.43185, 0.01311 


-0.86134, 0.58714 


3ct/ Ict 


5.13402, 0.00331 


3ct 


0.90169, 0.14221 


2ct 


CN 


Y 


(0.23,500) 


1.63419, 0.06861 


1.25677, 0.09989 


2ct/ 2ct 


3.68355, 0.01067 


3ct 


0.32701, 0.24697 


2ct 




Y 


(0.08,200) 


1.66239, 0.06672 


-1.02453, 0.63452 


2ct/ Ict 


-0.58847, 0.50465 


Ict 


1.26917, 0.09866 


2ct 


7 =1.3 




















CN 


Y 


(0.12,200) 


1.33772, 0.09214 


-0.33412, 0.42743 


2ct/ Ict 


-0.65832, 0.62597 


Ict 


0.1113, 0.3022 


2ct 


7 =1.4 




















CN 


Y 


(0.17,300) 


1.61382, 0.07001 


0.82904, 0.15277 


2ct/ 2ct 


1.65549, 0.06718 


2ct 


1.00266, 0.12868 


2ct 


7 =1.5 




















CN 


Y 


(0.23,200) 


1.73898, 0.06187 


0.60585, 0.18986 


2ct/ 2ct 


-0.82308, 0.57576 


Ict 


0.14452, 0.2937 


2ct 



